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Abstract 

Recent progress in understanding (2+l)-dimensional Yang-Mills {YM2+1) theory via the 
use of gauge-invariant variables is reviewed. Among other things, we discuss the vacuum 
wavefunction, an analytic calculation of the string tension and the propagator mass for 
gluons and its relation to the magnetic mass for YM^^i at nonzero temperature. (Talk 
given at the Workshop on Physical Variables in Gauge Theories, Dubna, Russia, September 
1999.) 
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1. Introduction 

In this talk I shall discuss a Hamiltonian approach to Yang-Mills theory in two spatial 
dimensions {YM2+1), where nonperturbative calculations can be carried out to the extent 
that results on mass gap and string tension can be compared with lattice simulations of the 
theory. The work I shall report on was developed over the last few years in collaboration 
with D. Karabali and Chanju Kim |]|, |2|, (Eventhough I shall concentrate on the pure 
YM2+1, some of the results can be extended to the Yang-Mills-Chern-Simons theory |^). 
Physical variables or gauge-invariant variables play a key role in our analysis, in keeping 
with the theme of this conference, and I hope this will be a nice example of how such 
variables can elucidate the nonperturbative structure of gauge theories. 

Before entering into the details of our work, let me say a few words about the relevance of 
YM2+1. Gauge theories (without matter) in (H-l) dimensions are rather trivial since there 
are no propagating degrees of freedom, although there may be global degrees of freedom 
on spaces of nontrivial topology. In (2-1-1) dimensions, gauge theories do have propagating 
degrees of freedom and being next in the order of complexity, it is possible that they provide 
a model simple enough to analyze mathematically and yet nontrivial enough to teach us 
some lessons about (3-|-l)-dimensional YM theories. Another very important reason to 
study YM2+1 is its relevance to magnetic screening in YM^-^-i at high temperature. Gauge 
theories at finite temperature have worse infrared problems than at zero temperature due 
to the divergent nature of the Bose distribution for low energy modes. A dynamically 
generated Debye-type screening mass will eliminate some of these, but we need a magnetic 
screening mass as well to have a perturbative expansion which is well defined in the infrared. 

A simple way to see that a magnetic mass can be dynamically generated is as follows. In 
the imaginary time formalism, with Matsubara frequencies ujn = 2TrnT, where T is the tem- 
perature, the gauge fields have a mode expansion as Ai{x, x^) = J2n ^i,n{x) exp{i2iTnTx^). 
At high temperatures and for modes of wavelength long compared to l/T, the modes 
with nonzero Matsubara frequencies are unimportant and the theory reduces to the the- 
ory of the = mode, viz., a three (Euclidean) dimensional Yang-Mills theory (or a 
(2-|-l)-dimensional theory in a Wick rotated version). Yang-Mills theories in three or (2-1-1) 
dimensions are expected to have a mass gap and this is effectively the magnetic mass of 
the (3-|-l)-dimensional theory at high temperature In order to incorporate this feature 
into YM^j^i at nonzero temperatures, one needs a decomposition of the YM^^i Feynman 
integrals wherein the YM2+1 pieces are isolated; in other words, one needs to identify the 
"slots" in the perturbative expansion of YM^^i where the YM2+1 results can be inserted. 
There is a recent analysis along these lines by Reinbach and Schulz 

Let me now start by recalling a couple of facts about YM2+i- The coupling constant 

has the dimension of mass and it does not run as the four-dimensional coupling does. 
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The dimensionless expansion parameter of the theory is k/e^ or e^//c, where A; is a typical 
momentum. Thus modes of low momenta must be treated nonperturbatively, while modes 
of high momenta can be treated perturbatively. There is no simple dimensionless expansion 
parameter. YM2+1 is perturbatively super-renormalizable, so the ultraviolet singularities 
are well under control. 

2. The parametrization of the fields 

Coming now to the details of our analysis, let us consider a gauge theory with group 
G = SU{N) in the ^0 = gauge. The gauge potential can be written as Ai = —if^Af, 
i = 1,2, where t° are hermitian N x A^-matrices which form a basis of the Lie algebra 
of SU{N) with [f^t^] = if'H'', i:T{tH^) = ^5"^. The spatial coordinates xi,X2 will be 
combined into the complex combinations z = xi —1x2, z = xi +1x2 with the corresponding 
components for the potential A = A-^ = \{Ai + iA2), A = Az = ^{Ai — iA2) = — (A^)^. 
The starting point of our analysis is a change of variables given by 

Az = -dzMM-^, Az = M^-^d-zM^ (1) 

Here M, are complex matrices in general, not unitary. If they are unitary, the potenial 
is a pure gauge. The parametrization (|^) is possible and is standard in many discussions of 
two-dimensional gauge fields. (There are also has similarities between (|l]) and the construc- 
tion of gauge-invariant particle states as discussed by McMullan, Lavelle and Horan at this 
workshop P].) A particular advantage of this parametrization is the way gauge transforma- 
tions are realized. A gauge transformation Ai — > A^ = Aig + dig , g{x) G SU{N) is 
obtained by the transformation M M^^^ = gM. The gauge-invariant degrees of freedom 
are parametrized by the hermitian matrix H = M^M. Physical state wavefunctions are 
functions of H. 

In making a change of variables in a Hamiltonian formalism, there are two things we 
must do: 1) evaluate the volume measure (or Jacobian of the transformation) which deter- 
mines the inner product of the wavefunctions and 2) rewrite the Hamiltonian as an operator 
involving the new variables. A consistency check would then be the self-adjointness of the 
Hamiltonian with the given inner product. We begin with the volume measure for the 
configuration space. 

3. The functional measure and inner product 

The YM Lagrangian in the Aq = gauge is given by 

L = J d^x 

where i?" = ^eij{diAj — djAf + f'^^'^A^Aj). By comparison of the kinetic term with the 
standard point-particle Lagrangian L = \g^vq^q^ , we see that the metric for the fields A, A 



e^ dAf dAf 1 

'2~dr~dr ~ 2^ 



(2) 
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is dsJ^ = J (fx 6Af5Af, with the corresponding volume dfi{A) = Y[x,a(^^^ix)dA"'ix), 
which is the standard Euchdean volume. From (|l|) we see that 

5A = -D{5MM-^) = - (d{5MM-'^) + [A,6MM~^]j 

5A = D{M^-^5M^) (3) 

which gives 

dfi{A) = {det DD) dn{M,M^) (4) 

where dfj,{M, M^) is the volume for the complex matrices M,M^, which is associated with 
the metric = 8 J Tt{6MM'^ M'^-'^SM^). This is given by the highest order differential 
form dV as dfi{M, Aft) = JJ^ dV{M, Aft) where 

dV{M,M^) oc ea,...aSdMM'^)a, A ... A {dMM-^)a^ 

X ebi...b„(Mt-idMt),^ A ... A {M^-^dM%^ (5) 

where n = dimG = dmiSU {N) = N'^ — 1. (There are some constant numerical factors which 
are irrelevant for our discussion.) The complex matrix M can be written as M = Up, where 
U is unitary and p is hermitian. This is the matrix analogue of the modulus and phase 
decomposition for a complex number. Since gauge transformations act as M ^ M^^^ = gM, 
we see that U represents the gauge degrees of freedom and p represents the gauge-invariant 
degrees of freedom on M. Substituting M = Up, (P) becomes 

dV{M,M^) (X ea,...aAdpp'^ + P'^dp)a, A ... A {dpp-^ + p-^dp)a,, 
xeb,...bAU-^dU)b,A...A{U-^dU)b„ 
oc ea,...aSH-^dH)a,A...A{H-^dH)aJp{U) (6) 

Here dp{U) is the standard group volume measure (the Haar measure) for SU{N). Upon 
taking the product over all points, dp[U) gives the volume of the entire gauge group (namely 
all S'J7(A^)-valued functions) which we denote by vol{Q^) and thus 

dp{M, Aft) = W dV{M, Aft) vol{g^) = dp{H) vol{g^) (7) 

X 

dn{H) = Y[ det r [dip''] 

x,a 

detrH^V'" =ea,...aAH-^dH)a,...{H-^dH)a^ (8) 

a 

We have parametrized H in terms of the real parameters and H^^dH = dip^rak{^)tk- 
The volume element or the integration measure for the gauge-invariant configurations can 
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now be written as 



vol{G*) 



vol{G*) 
(det D,D,) 



dn{M, Mt) 



{det DD)dfi{H) 



(9) 



where we have used (^) . The problem is thus reduced to the calculation of the determinant 
of the two-dimensional operator DD. This is well known The simplest way to evaluate 
this is to define T = log det DD, which gives 

6r 



-i Tr 



(10) 



{T^)mn = —ifmn ^'^^ the generators of the Lie algebra in the adjoint representation. The 
coincident-point limit of D~^(x,y) is singular and needs regularization. With a gauge- 
invariant regulator, one finds 



Tr 



2CA 



Tr 



{A - M^-^dM^)e 



{ID 



where ca^"'' = p-^^-fmn. gq^^j ^ fo^, SU{N). Using this result in (jl^) and 

integrating we get 



{det DD) 



det' dd 



n dimG 



exp [2cA SiH)] 



(12) 



S{H) is the Wess-Zumino-Witten (WZW) action for the hermitian matrix field H given by 



1 



27r 



127r 



S{H) = — / TiidHdH-') + / €'"'''Ti{H-'d^HH-'duHH-'daH) (13) 



We can now write the inner product for states |1) and |2), represented by the wavefunc- 

(14) 



tions ^'i and ^'2, as |]To| 



(1|2) = J dii{H)e^^^ ^^"^ ^\^2 

4. Transforming the Hamiltonian 

The next step is the change of variables in the Hamiltonian. However, there is some 
further simplification we can do before taking up the Hamiltonian. We would expect the 
wavefunctions to be functionals of the matrix field H, but actually we can take them to be 
functionals of the current of the WZW model ( p^ ) given by J = {cA/T^)dzH H^^. First of 
all we notice that the Wilson loop operator can be constructed from J alone as 



WiC) = TrP e"^c 



(Adz+Adz) 



(15) 
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Since the Wilson loop operator can provide a complete description of gauge-invariant ob- 
servables, it is sufficient to take wavefunctions to be functions of J. There is also a conformal 
theory argument to show that it is sufficient to conside only functions of J ||l|, • 

This means that we can transform the Hamiltonian Ti. = T + V to express it in terms 
of J and functional derivatives with respect to J. This is achieved by the chain rule of 
differentiation 



— I E?E? ^ 



V 



1 

2^2 



6J^{u) 



5^' 



6Af{x)6AI{x) SJ^iu] 



+ 



6J^{u) 6J^{v) 



<5J"(ii) 6J\v 



(16) 



Regularization is important in calculating the coefficients of the two terms in T. Carrying 
this out we find 



m 



V 



TT 



mcA 



dJa{x)dJa{x) 



5J<'{u) 5J^{v) 



where m = c^ca/'^t^ and 



rLab{u,v) 



CA 



JabcJ^iv) 



(17) 

(18) 
(19) 



vr^ (n — v)'^ tt{u — v) 

The first term in T shows that every power of J in the wavefunction gives a value m to 
the energy, suggesting the existence of a mass gap. The calculation of this term involves 
exactly the same quantity as in ( p^ ) and with the same regulator leads to (p!7[), i.e., 

-2 " S'Jaix) 



d'y- 



(?CA 



Aft Tr 

27r 



TO-^(y,x 



8A\y)bA\y) 

= m Ja{x) (20) 

Finally, (17,1^) (with regularizations taken account of) give a self-adjoint Hamiltonian 
which, as I mentioned before, is a nice consistency check. 

5. An intuitive argument 

The next step is to solve the Schrodinger equation, at least for the vacuum state and 
the low lying excited states. However, before taking this up, I shall give a short intuitive 
argument for the existence of a mass gap. 
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First of all notice that the total volume of the configuration space as given by (^) is 
finite, modulo regularization of the Laplacian dd, and can be written as fin, ffl] 



vol{Q^) 



jcPx 



n -dimG 



(21) 



In contrast to this, the corresponding result for an Abelian theory (which has ca = 0) 
is infinite. (If we make a mode decomposition of H or over the eigenmodes of the 
Laplacian dd, the integration over the amplitude of each mode is finite for the nonabelian 
case because of the exponential; the divergence arises from the infinity of modes and can be 
regularized by truncation to a finite number of modes. For the Abelian case, the integration 
for each mode is divergent.) This result is encouraging as regards the question of the mass 
gap. One can go further and make a slightly better argument. The crucial ingredient is 
the measure of integration in the inner product (|l^). Writing A£^, AS for the root mean 
square fluctuations of the electric field E and the magnetic field B, we have, from the 
canonical commutation rules = —idijd""^, AE AB ~ k, where k is the momentum 

variable. This gives an estimate for the energy 



^=ohri^ + -^ (22) 




For low lying states, we minimize £ with respect to AB^, AB^^^ ~ e^k, giving £ ^ k. This 
corresponds to the standard photon. For the nonabelian theory, this is inadequate since 
{H) involves the factor e^'^-^ ^ ^^^^^ 

{H) ~ j dfi{H)e^^^^^^^ ^{e^E^ + B'^/e^) (23) 

In terms of B, the WZW action goes like S{H) [— (c^/vr)^ / B{l/k'^)B + ...]; we thus see 
that B follows a Gaussian distribution of width AB'^ ~ irk"^ /ca, for small values of k. This 
Gaussian dominates near small k giving AB'^ ~ A;^(7r/c^). In other words, eventhough £ 
is minimized around AB^ ~ k, probability is concentrated around AB^ ~ k'^iir/cA)- For 
the expectation value of the energy, we then find £ ~ c^ca/'^t^ + 0{k'^). Thus the kinetic 
term in combination with the measure factor e^'^-*'^*^^^ could lead to a mass gap of order 
e^CA- The argument is not rigorous, but captures the essential physics as we shall see in a 
moment. 

6. The vacuum wavef unction 

Let us now consider the eigenstates of the theory. The vacuum wavefunction is pre- 
sumably the simplest to calculate. Ignoring the potential term V for the moment, since 
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T involves derivatives, we see immediately that the ground state wavefunction for T is 
<I>o = 1- This may seem like a trivial statement, but the key point is that it is normalizable 



with the inner product (|T4|); in fact, the normalization integral is (21). Starting with this, 
we can solve the Schroodinger equation taking \I'o to be of the form exp(P), where P is a 
perturbative series in the potential term V (equivalent to a 1/m-expansion). We then get 



P 



TT 



m?CA 



Tr / : dJdJ 



^— ) Tr / [: dJ{Vd)dJ + \dJ[J,d\j] :] 

^— ^Tr f[:dj{vdfdj + -[vdj, dj]d'^j + -[vd^j, j]d^j 

m^CA/ 7 9 9 

- ^[J, dj][dj, d^j] - ^[j,dJUd'j] -HOi^) 



(24) 



where Dh = (c^/7r)(9/i — [J, h]. The series is naturally grouped as terms with 2 J's, terms 
with 3 J's, etc. These terms can be summed up; for the 2J-terms we find 



= exp 



1 

2^ 



Ba{x) 



1 



(m + \/ 



V2) 



Ba{y) +0(3J) 



J x,y 



(25) 



The first term in (25) has the correct (perturbative) high momentum limit, viz.. 



^'o ~ exp 



Jx.y 



1 



Ba{y) +0{3J) 



x,y 



(26) 



Thus although we started with the high m (or low momentum) limit, the result (|2^) does 
match onto the perturbative limit. The higher terms are also small for the low momentum 
limit. 

We can now use this result to calculate the expectation value of the Wilson loop operator; 
for the fundamental representation, it is given by 



{Wf{C)) 



constant exp [— cr^c*] 

^ im - 1 



Svr 



(27) 



where Ac is the area of the loop C . a is the string tension. This is a prediction of 
our analysis starting from first principles with no adjustable parameters. Notice that the 
dependence on and is in agreement with large-A^ expectations, with o depending 
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only on the combination e N as N — > oo. (The first correction to the large- limit 
is negative, viz., — {e^N)/2N'^^/^ which may be interesting in the context of large-A^ 
analyses.) Formula (|^ gives the values ^/a/e^ = 0.345, 0.564, 0.772, 0.977 for iV = 2, 3, 4, 5. 
There are estimates for a based on Monte Carlo simulations of lattice gauge theory. The 
most recent results for the gauge groups SU{2), SU{3), 5C/(4) and SU{5) are \fole^ = 
0.335, 0.553, 0.758, 0.966 We see that our result agrees with the lattice result to within 
~ 3%. 

One might wonder at this stage why the result is so good when we have not included 
the 3J- and higher terms in the the wavefunction. This is basically because the string 
tension is determined by large area loops and for these, it is the long distance part of the 
wavefunction which contributes significantly. In this limit, the 3 J- and higher terms in (|25| ) 
are small compared to the quadratic term. 

We have summed up the 3J-terms as well. Generally, one finds that P, when expressed 
in terms of the magnetic field, is nonlocal even in a (l/m)-expansion, contrary to what one 
might expect for a theory with a mass. This is essentially due to gauge invariance combined 
with our choice of = Oj it has recently been shown that a similar result holds for the 
Schwinger model [12|. 

7. Magnetic mass 



I shall now briefly return to the magnetic mass. From the expression (17) we see 
immediately that for a wavefunction which is just J^, we have the exact result T J"" = m J"". 
When the potential term is added, J" is no longer an exact eigenstate; we find 



(T + V) r = Vm^ - V2 r + • • • (28) 

showing how the mass value is corrected to the relativistic dispersion relation. 

Now J" may be considered as the gauge-invariant definition of the gluon. This result 
thus suggests a dynamical propagator mass m = 6^0^/271 for the gluon. A different way 
to see this result is as follows. We can expand the matrix field J in powers of ipa which 
parametrizes H, so that J ~ {cA/T^)dipata- This is like a perturbation expansion, but a 
resummed or improved version of it, where we expand the WZW action in exp(2cyi5(ff)) 
but not expand the exponential itself. The Hamiltonian can then be simplified as 



^ ^ H [-Z^ + - V^)<Pa{x)] + ... (29) 



where (j)a{k) = \J CAkk/{2'Km)i^a{^), in momentum space. In arriving at this expression we 
have expanded the currents and also absorbed the WZW-action part of the measure into 
the definition of the wavefunctions, i.e., the operator (|29|) acts on ^ = e'^^'^'^^^^'. The above 
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equation shows that the propagating particles in the perturbative regime, where the power 
series expansion of the current is appropriate, have a mass m = e^CA/27r. This value can 
therefore be identified as the magnetic mass of the gluons as given by this nonperturbative 
analysis. 

For SU{2) our result is m ~ 0.32e^. There have been two lattice calculations of the 



propagator mass of gluons for SU{2); the values are m k, 0.35e^ and 0.46e^ jl^, 14|. There 
is reasonable evidence for a 'constituent gluon' picture for glueballs in YM^ from lattice 
analysis, so another suggestion has been to extract a constituent mass for a gluon and 
interpret it as the magnetic mass This gives a value m ~ 0.31e^ — 0.40e^. Considering 
the difficulties of a lattice estimate and the variance within these calculations, we cannot 
draw any definite conclusion in comparing with our analysis. 

8. Excited states 

Eventhough J is useful as a description of the gluon, it is not a physical state. This 
is because of an ambiguity in our parametrization (|l|). Notice that the matrices M and 
MV{z) both give the same A, A, where V{z) only depends on z and not z. Since we have 
the same potentials, physical results must be insensitive to this redundancy in the choice 
of M; in other words, physical wavefunctions must be invariant under M MV{z). J is 
not invariant; we need at least two J's to form an invariant combination. An example is 



*2=/ f{x,y)[dJa{x){H{x,y)H'\y,y))abdJb{y)\ (30) 

Jx,y 

This is not an eigenstate of the Hamiltonian. Since there are two J's we should expect 
at least a mass of 2m, but beyond that it is difficult to say anything very conclusive, see 
however [^. 

This work was supported in part by a grant from the National Science Foundation. 
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